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Method and system for decomposition of multiple channel signals 

5 

FIELD OF THE INVENTION 

This invention relates to instrumentation and related methodology for 
measurements of an optic composite noisy signal, and in particular, for tracking and 
recovering multiple sinusoids from a multiple-wavelength signal generated by a 
10 photodetector. 

BACKGROUND OF THE INVENTION 

The problem of tracking and recovering multiple sinusoids from an optic 
multiple-wavelength signal is of high practical importance. For example, this 
problem arises in the non-invasive measurements of various blood parameters, e.g., 

15 blood oxygenation, that employ the detection of light transmitted or reflected from 
flie location on the patient's body under measurement, and indication of the 
presence of various blood constituents based on known spectral behaviors of these 
constituents. Most of these techniques utilize a measurement optical device or 
probe attached to the patient's finger, which includes an optical assembly including 

20 a block of light emitting devices (e.g., LEDs) for irradiating the finger with light, 
and a photodetector (PD) for detecting the finger's light response. One of the 
conventional devices of the kind specified, such as a pulse oximeter, which is the 
generally accepted standard of everyday clinical practice, provides for measuring 
enhanced optical pulsatile signals caused by the changes Jn the volume of a blood 

25 flowing through a fleshy medium (e.g., finger). In conventional pulse oximeters as 
well as in other plediysmographic devices, an optical, amplitude-modulated light 
passes through die investigated sample (e.g., finger) at several (usually 2-3) 
different carrier wavelength (channels). Thereafter, a relationship between the 



i 

-2- 

intensities of the corresponding light responses is measured for determining the 
corresponding blood parameters. This technique utilizes the fact that different 
wavelength optic signals (channels) are transmitted differently, when passing 
through Ihe sample, in accordance with the scattering and/or absorbing properties 
5 of the blood. In such techniques, the intensity of the signal obtained from the blood 
measurements (e.g., the intensity of photo-plethys olographic signal (PPS)) for the 
particular optical wavelength is obtained by means of a modulation-demodulation 
(M-D) procedure. For this purpose, the light emitting devices are switched on and 
off at a presumed rate (e.g., about 3-10 kHz), providing thereby a sequence (train) 
io of light pulses for modulation of the corresponding carrier waves. The total 
pulsatile optical signal is detected by the photodetector and demodulated. The 
wavelength-related intensities are then derived from the demodulated optical 
signals. 

In such conventional techniques, LEDs usually operate in one of two 
15 regimes, serial or parallel. In the serial mode regime, the operation of LEDs is 
separated in time. On the other hand, in the parallel mode regime, the LEDs are 
activated concurrently, but with the diverse modulation rates, or, in other words, the 
LED signals are separated in the frequency domain. The demodulation procedure 
used for separating different optic channels from each other is usually based on 
20 employing a bank of frequency-selective filters arranged either in cascade or 
parallel configurations. 

More specifically, in the serial operation mode regime, each LED oscillates 
over a certain, relatively long period (including about 10-30 on-off LED ignition 
cycles), while the other LEDs are kept inactive over this period An output signal, 
25 sensed by the photodetector, is then demodulated and associated with the 
corresponding wavelengths. Such a procedure is repeated for the next LED, et 
cetera. 

In fact, the complexity of the modulation-demodulation technique rapidly 
increases with the number of the desired optical channels. Because of this, a 
30 majority of the commercial oximeters utilize few, e.g., 2-3 wavelengths. When such 
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a relatively small number of optical channels is used, one can conduct PPS 
measurements for different optic channels in series, one by one. In such a case, as 
described above, each LED is activated serially over a certain period, used for 
demodulation, thus outputting a modulated signal related to the particular 
5 wavelength. 

However, there are several critical points behind the serial mode regime. As 
the number of channels increases (e.g., close to ten), then the serial order comes 
into conflict with the required access rate. In other words, the physical and 
chemical processes, e.g., aggregation of red blood cells, occurring in the measured 
10 sample can be different for the time periods corresponding to the activation of 
different LBDs. 

One of the shortcomings associated with the serial operation mode regime is 
the phase shift between the optical channels. Thus, the trai n of pulses (provided by 
a generator for LEDs* activation) switches between different LEDs by means of a 
15 multiplexer (MUX) such that only a single LED operates at a time. To this end, the 
longer LED operation period the higher delay between the different channels. 

On the other hand, decreasing periods of LED activity restricts the 
demodulator selectivity and accuracy. So, in the serial mode, the demodulator 
selectivity and accuracy are in contradiction with the measuring process 
20 simultaneity. 

Moreover, the serial mode regime yields various problems in signal 
processing. In this connection, it should be noted that the signal obtained from the 
blood measurements is usually affected by the electrical 50 Hz interference, 
lamp-induced 100 Hz optical noise, RE, and other disturbances. Hence, in the serial 

25 operation mode, the continuous interference and trends, that affect the output of the 
photodetcctor, cannot be rejected properly, because each channel, containing a part 
of the interference, is processed periodically and independently. A conventional 
prefihering of the photodetector's output (in order to enhance the integral signal 
before the channel separation) breaks the independence between the consequent 

30 signal frames. 



I 

-4- 

Replacing the serial order in optic by a parallel order in the conventional 
techniques leads to classical frequency-separation methods. In the parallel 
operation mode regime, all LEDs arc activated concurrently, however with different 
switching rates. The further recovering of different optic channels from the 
5 composite signal (when the number of channels is small) is usually performed by 
means of a frequency-selective filter-bank. 

It should be noted that in parallel mode regime, prefiltering as well as other 
conventional signal enhancing methods may generally become feasible. 
Nevertheless, in spite of the flexibility and universality, this regime becomes 
10 impractical when the number of wavelengths increases significantly, e.g. more than 
three. 

For example, in the parallel operation mode regime the frequency spectrum 
lines (optic channels) should preferably be placed as far as possible from each 
other. This requirement can be hardly satisfied in practice for the case of a large 

15 number of LEDs. More specifically, a modulation-demodulation concept should 
meet the given technical restrictions. In practice, the modulation process is formed 
by triggering the LEDs at a certain rate. Thus, on the one hand, the LED triggering 
rate is restricted, for example, by the transient of the photodetector circuitry. On the 
other hand, it should be taken into account that the resulting signal may be of a 

20 rather complex form and a care should be taken that the spectra do not intersect. In 
particular, since the output LED signals are non-sinusoidal (i.e., nearly square 
waves), some harmonics of the low-rate LEDs may interfere with the fundamentals 
belonging to high-rate LEDs. In other words, the first and second harmonics of the 
slowest L ED d ctcrm i ne t he i nterval o f p ossible m odulation r ates, H ence, a 11 the 

25 other triggering rates should be distributed properly within this restricted interval. 
For instance, i f the lowest modulation rate is defined at 2000 Hz, then (for example 
for the case often optic channels) the remaining nine on-off LEDs' rates should be 
placed in the interval between 2000 Hz and its second harmonic 4000 Hz. Thus, for 
the nine channels the frequencies can be placed over equally spaced 200 Hz 

30 increments at points such as: 2200, 2400, 3800 Hz. It should therefore be 
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apprcciated that as a number of the optic channels increases, the frequency distance 
between the modulated signals decreases, thereby causing the interference problem. 
Thus, in the case of multiple optic channels, the parallel operation mode may result 
in a "crowd" signal with closely spaced and strongly coupled components. 
5 Moreover, due to the deviations from linearity in the signal acquisition and 

transformation tract all harmonics may interfere in the multi-channel system, thus 
producing the difference and sum of the frequencies. In practice, the signal 
spectrum comprises not only the first and higher harmonics, but their combinations 
as well, i.e., multiple sub-harmonics. The latter may fall close to the fundamental 
10 frequency, thus increasing the complexity of the signal decomposition in the 
parallel mode. 

Next, the parallel mode regime may yield a corresponding extra complexity 
of the circuitry of the multi-channel system, because the number of the wires 
coupled to the probe grows in direct proportion with the number of the optic 
15 channels. 

Moreover, the c oncurrent r un o f s everal L EDs c an c ause t he o verloading 
(saturation) of the photodetector, thus restricting the signal-to-noise ratio for each 
particular wavelength. 

Another difficulty associated with utilization of the conventional parallel 
20 mode regime in signal decomposition techniques is in the fact that the sensed 
optical intensity may crucially differ for different channels. Thus, the channel 
separation problem can be compounded by the interference between the closely 
spaced strong and weak signals. 

Due to the aforementioned disadvantages and also other known in the art 
25 reasons, a utilization of the pure parallel mode regime in multi-channel optical 
systems is impractical. 

Concluding the above consideration, neither pure serial operation mode 
regime nor parallel mode regime alone is a reasonable and practical candidate for 
the system that intend to use a large number of wavelengths. Hence, when a 



I 



~6- 

multi-channel signal acquisition is of interest, a more sophisticated technique 
should be advantageous. 



SUMMARY OF THE INVENTION 

5 Thus, there is still a need in the art for, and it would be useful to have a 

novel method and a system for tracking and recovering multiple sinusoids from a 
multi-channel signal reflecting characteristics of a blood perfused fleshy medium. It 
would be advantageous that such a technique could decompose the signal having a 
high level of optical disturbances and trends. Moreover, the technique should, inter 
10 take into account the fact that the output waveforms of the light sources (e.g., 

LEDs) are not ideally square and may vary in shape from LED to LED. Since these 
shapes are not fixed and may vary due to the light intensity, temperature, etc., the 
technique should take into account that the signal sensed by a photodetector can be 
rich in harmonics that cannot be ignored. 
15 It would also be desirable, that the system for the measurements of blood 

parameters be properly adjusted, in order to combine several requirements, such as, 
high modulation rate, continuity of measurements for each wavelength, 
simultaneity, selectivity, etc. 

The present invention satisfies the aforementioned need by providing a 
20 novel measuring system for determining blood parameters of a blood perfused 
fleshy medium. The system includes a generator, a multiplexer (MUX) coupled to 
the generator, a measuring probe coupled to the MUX and an analyzer coupled to 
the probe and, preferably, to the MUX. The generator is for providing a train of 
activating pulses. The MUX is configured for switching said activating pulses 
25 between different optic channels. The MUX is adapted for providing a composite 
mode regime for the activating pulses. The composite mode regime represents a 
combination of parallel and short serial modes of the activating pulses. The probe 
includes an illumination assembly having a plurality of light sources coupled to the 
MUX and activated by the activating pulses that provide ignition of the light 



sources. The light sources are configured for generating amplitude-modulated light 
of more than two different optic channels having wavelength in a range where the 
scattering properties of blood are sensitive to light radiation. The probe also 
includes a photodetcctor adapted for sensing a light response of the medium and 
5 generating a multiple channel signal reflecting blood characteristics. 

The analyzer is configured for analyzing the multiple channel signal 
generated by the photodetector. The analyzer includes an analog-to-digital 
converter, a first decimator, a digital signal processor having an adaptive resonator 
bank unit, an output filtering unit, and a second decimator unit. The adaptive 
10 resonator bank unit is configured for (i) filtering the multiple channel signal out of 
trend signals, optical and electrical disturbances along with noise and (ii) separating 
the multiple channels from each other. The adaptive resonator bank unit is a 
closed-loop resonator bank, based on the spectral observer scheme, and coupled to 
the frequency adaptation loop. The output filtering unit is configured for obtaining 
15 synchronous in time evolutions of the light responses of the medium for the 
multiple channels. The second decimator is designed for outputting the time 
evolutions of the light responses at a lower sampling rate (e.g., 25-30 sample/sec), 
which may be required for certain applications. 

The present invention also satisfies the aforementioned need by providing a 
20 method for decomposition of a multiple channel signal reflecting characteristics of 
a blood perfused fleshy medium. The method can, inter alia, be used for 
determination of desired blood parameters. 

The method includes the step of illuminating a portion of the medium by 
amplitude-modulated light of more than two different optic channels. The light has 
25 wavelengths in a range where the s cattering properties of blood are sensitive to 
light radiation. The amplitude-modulated light is activated in a composite mode 
regime representing a combination of parallel and serial modes of ignition of the 
light sources. 
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According to one example, the composite mode regime is a short serial 
mode regime, where each light source is modulated over a short interval, e.g., 1-5 
on-off light ignition cycles. 

According to another example, the composite mode regime is a short 
5 serial-parallel mode regime. In this regime, the light sources can be divided into 
two or more groups, where each group is operated in the short serial mode, 
however, with a different rate. 

According to yet another example, the composite mode regime is a 
mixed-rate short serial mode regime. In this regime, at least two pulse trains 
10 having different rates are initially mixed, thereby resulting in a single 
complex-frame signal, and then applied to all the light sources at a certain time 
moment in the short serial mode regime. 

The method also includes the step of sensing a light response of the medium 
and generating the multiple channel signal reflecting blood characteristics. Further, 
15 the method includes the step of analyzing the multiple channel signal- In turn, the 
analyzing of the multi-channel signal includes analog-to-digital conversion and 
decimating the signal. Moreover the analyzing of the multiple channel signal 
includes filtering the multiple channel signal and separating the multiple channels 
from each other. Both the filtering and separating are carried out by applying the 
20 adaptive resonator bank to the multiple channel signal The analyzing also includes 
providing time evolutions of the 1 ight responses of the medium for the multiple 
channels. 

Finally, the method comprises the step of deriving at least one blood 
characteristic parameter as a relation between the different time evolutions of the 
25 light responses (light intensity signals) of the medium obtained by filtering sinusoid 
signals corresponding to the optic channels. 

There has thus been outlined, rather broadly, the more important features of 
the invention in order that the detailed description thereof that follows hereinafter 
may be better understood. Additional details and advantages of the invention will 
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be set forth in the detailed description, and in part will be appreciated from the 
description, or may be learned by practice of the invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 

5 In order to understand the invention and to see how it may be carried out in 

practice, a preferred embodiment will now be described, by way of non-limiting 
example only, with reference to the accompanying drawings, in which: 

Fig. 1 is a schematic block diagram of a measuring system for the 
determination of blood parameters is illustrated, according to one embodiment of 
io the invention; 

Fig. 2 is a schematic block diagram illustrating the closed-loop resonator 
bank configuration; 

Fig. 3 is a schematic block diagram illustrating a synchronization between 
the light sources and resonators in the system of the present invention; 
is Fig. 4 is a schematic block diagram of a system of the present invention 

operating in a short serial mode regime, according to one example of the invention; 

Figs. 5 illustrates an example of wavefoims produced by different elements 
of the system shown in Fig, 4; 

Fig. 6 is a schematic block diagram of a system of the present invention 
20 operating in a short serial-parallel mode regime, according to another example of 
the invention; 

Figs. 7 illustrates an example of waveforms produced by different elements 
of the system shown in Fig. 6; 

Fig. 8 is a schematic block diagram of a system of the present invention 
25 operating in a mixed-rate short serial mode regime, according to yet another 
example of the invention; and 

Figs. 9 illustrates an example of waveforms produced by different elements 
of the system shown in Fig. 8, 
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DETAILED DESCRIPTION OF THE INVENTION 

The principles and operation of the process and system according to the 
present invention may be better understood with reference to the drawings and the 
accompanying description, wherein like reference numerals have been used 
5 throughout to designate identical elements. It is understood that these drawings are 
given for illustrative purposes only and are not meant to be limiting. It should be 
noted that the blocks in the drawings illustrating various embodiments of the 
system of the present invention are intended as functional entities only, such that 
the functional relationships between the entities are shown, rather than any physical 
i o connections and/or physical relationships. 

Some portions of the detailed descriptions, which follow hereinbelow, are 
presented in terms of algorithms and symbolic representations of operations on data 
represented as physical quantities within registers and memories of a computer 
system. An cilgorithm is here conceived to be a sequence of steps requiring physical 
15 manipulations of physical quantities and leading to a desired result. Usually, 
although not necessarily, these quantities take the form of electrical or magnetic 
signals capable of being stored, transferred, combined, compared, and otherwise 
manipulated In the present description, these signals will be referred to as values, 
elements, symbols, terms, numbers, or the like. Unless specifically stated otherwise, 
20 throughout the description, utilizing terms such as "processing* 1 or "computing" or 
"calculating" or "determining" or the like, refer to the action and processes of a 
computer System, or similar electronic computing device, that manipulates and 
transforms data. 

Referring to Fig. 1, a schematic block diagram of a measuring system 10 for 
25 determination of blood parameters is illustrated, according to one embodiment of 
the invention. The system includes a generator ll^a multiplexer (MUX) 12, a probe 
13 and an analyzer 14 coupled to the probe 13. The generator 11 provides a series 
of pulses. These pulses axe fed to the MUX 12, which is a switch configured to 
provide a composite operation mode regime of activation the probe 13. 
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Thc probe 13 comprises an illumination assembly 131 activated by the 
pulses generated by the generator and designed for generating light of different 
wavelengths (more than two different wavelengths). The probe 13 further includes 
a photodetector 132 adapted for sensing a light response of a blood perfused fleshy 
5 medium 130, and generating a multiple channel signal reflecting characteristics of 
the medium. Examples of the blood perfused fleshy medium include, but are not 
limited to, a patient's finger, a toe, a wrist, an elbow, etc. The illumination assembly 
131 includes a plurality of light sources (e.g., LEDs), each associated with a 
wavelength selected in a range where the scattering properties of blood are 
10 sensitive to light radiation. When required, the illumination assembly 131 can be 
driven by a controllable drive mechanism that is not specifically shown. Examples 
of the photodetector 132 include, but are not limited to, spectrophotometers and/or 
photodiodes typically equipped with amplifying means, which are not specifically 
shown. 

15 According to the invention, the MUX 12 is configured to provide a 

composite mode regime for ignition of the light sources. In other words, the light 
sources of the probe 13 can be run in the mode regime being different from the 
conventional pure parallel or serial mode regimes (described in the background 
section of the application). According to the invention, a mode regime representing 

20 a combination of parallel and short serial modes is employed, where pulses may 
have different rates. According to one embodiment of the invention, each LED is 
modulated over a short interval, e.g., over 1 to 5 on-off ignition cycles. This mode 
hereinafter will be referred to as a short serial mode regime (in contrast to the 
conventional serial mode regime in which each LED alternates over a considerably 

25 long interval, usually more than 10 on-off ignition cycles). 

According to another embodiment of the invention, the LEDs can be 
divided into two or more groups, where each group is operated in the short serial 
mode, however, with a different rate. This mode regime hereinafter will be referred 
to as a short serial-parallel mode regime. According to this embodiment, the 
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number of LEDs operating either in parallel or in series may be controliably 
adjusted. The duration of the LED's activity becomes a design parameter as well. 

According to yet another embodiment of the invention, two or more pulse 
trains provided by the generator are initially mixed, thereby providing a single 
5 complex-irame signal, which is further applied to all the LEDs in a short 
serial-parallel mode regime. This mode regime hereinafter will be referred to as a 
mixed-rate short serial mode regime. According to this embodiment, each LED 
switches in a short serial mode by a triggering signal that represents a mixture of at 
least two generator waveforms of different rates- 
10 The optical measurements are carried out to detect the light response of the 

medium (light transmitted through the medium or the light reflected therefrom, as 
the case may be), and to generate data representative thereof (constituting sensed 
signal). 

The photodetector 132 of the probe 13 is coupled to the analyzer 14 
15 responsive to the sensed signal. The multiple channel signal generated by (he 
photodetector 132 is a composite signal including multiple sinusoids (having 
time-varying amplitudes) related to the LEDs' activity along with signal trends, 
disturbances and noise. The analyzer 14 is a computer device having an 
analog-to-digital converter (ADC) 140 5 a digital signal processor (DSP) 141 and a 
20 monitor (not shown) for presenting the measurement results. This analyzer 14 may 
be specially constructed for the required purposes, or it may comprise a general 
purpose computer selectively activated or reconfigured by a software program 
stored in computer's memory (not shown), such as, but is not limited to, any type of 
disk including floppy disks, optical disks, CD-ROMs, and magnetic-optical disks, 
25 read-only memories (ROMs), random access memories (RAMs), EPROMs, 
EEPROMs, magnetic or optical cards or any type of media suitable^ for storing 
electronic instructions, and each coupled to a computer system bus. The DSP 141 
includes, a first decimator unit 142, an adaptive resonator bank unit 143, an output 
filter unit 144 and, optionally, a second decimator 145. 
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The ADC unit 140 is a known unit designed for sampling and digitizing the 
signal obtained from the photodetector 132. The first decimator 142 is configured 
for reducing the sampling rate of the digital signal treated in the DSR The various 
examples of the decimating process utilized by the first decimator will be described 
5 below. The second decimator is designed for outputting the time evolutions of the 
light responses at a lower sampling rate (e.g., 25-30 sample/sec), which may be 
required for certain applications. 

It should be appreciated that the DSP 141 is preprogrammed by a 
mathematical model capable of analyzing the received output (sensed signal) of the 
10 detection assembly and determining desired parameters of the patient's blood. 

The adaptive resonator-bank (ARB) 143 constitutes of frequency-selective 
filters - resonators, configured for tracking and recovering the multiple 
time-varying sinusoids from the composite noisy signal generated by the 
photodetector 132, The ARB 143 concurrently performs the amplitude 
15 demodulation for each channel along with frequency tracking and_rejecting of the 
undesirable spectral lines, such as trends, disturbances and noise. 

Finally, as a result of the operation of the output filtering unit 144, the light 
intensity signals associated with the LEDs' selected wavelength can be obtained 
and used for deteimining blood parameters by using an appropriate known method 
20 based on a relationship between these light intensities. For example, methods 
described in U.S. Pat. Nos. 6,400,972 and 6,587,704 assigned to the Applicant of 
this application, and incorporated herein by reference can be used for the 
determining of blood parameters. 

According to the invention, the adaptive resonator bank (ARB) is based on 
25 the closed-loop resonator bank representing a parallel connection of narrow-band 
resonance filters with a feedback loop. The general mathematical principles oFthe 
application of the adaptive resonator bank (ARB) to the composite signal generated 
by the photodetector 132 will be expounded hereinbelow. 

Generally, the composite signal generated by the photodetector 132 can be 
30 represented by a time-varying vector s=s(/*)=(ai,.-., <ryv) T , where ac=Apos(cQin+ <p$ 
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is the z-th sinusoid having the amplitude A h the angular frequency co L and the phase 
(pi\ n is the discrete time. For example, the composite signal can be represented by a 
sum of sinusoids corrupted by the noise v—v(n) 9 to wit, y=ai+. « .+cnv +v. A retrieval 
of the sinusoids is usually performed by applying a filter bank that generates an 

5 M<1 output vector y={y\ y^), where each component j// (i=l,...,N) relates to the 

corresponding input element a/. 

According to the invention, the determining of the vector y={y\ yv) is 

carried out in terms of the state-space formulation of the Kalman-like spectral 
observer (see, for example, GC Goodwin, et aL, "Sinusoidal disturbance rejection 
10 with application to helicopter flight data estimation" IEEE Trans. Acoustics, 
Speech, and Signal Processing, 1986, V 34(3), PP. 479-484), which for a single 
tone coi can be presented by the following closed-loop pattern: 

*<*)•■ >>(*)-*,*! 00 

where x?=xfa)=[cj si] is the 2x1 state vector comprising the in-phase and 
15 quadrature components ci and st, respectively; G, is the 2x1 gain vector; hi =[1 0] is 
the observation matrix, and 7} is the 2x2 transition matrix 

r '{-% 3 < 2 > 

where Q=cos(co,) and S i =sin(a,,) ) f=l, ...,N. 

Accordingly, the AT-sinusoid Spectral Observer can be represented by 
x(« + 1) m Tx(n) + Ge(«), 

where *-j{iiM*i.—, x,v] T ; G=lG it ..., Gwf; h=[h\,... y h N \ and the 2Nx2N 
transition matrix is defined as T=blockdiag{Ti} y i=l. .... N. 

The transfer function (TF) connecting the prediction error e with the 
particular /-th output x h F,{z):e^xi , corresponding to the z'-th sinusoid can be 
25 obtained by 
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where z is the forward time-shift operator. 

For example, in the case when the gain is represented by Gr z ygr='j{Cu -iS/J , 
where y is the relaxation factor, the corresponding mathematical model yields the 
following expression for the /-th sinusoid TF: 

y [CgZ-1 

5 F = 

1 z 2 -2C,-z-fl|_ 

As can be appreciated, the numerator of the fraction in Eq. (5) matches the 
corresponding sinusoidal oscillator. In other words, each particular component of 
the processed optical signal is associated with the corresponding resonator in the 
ARB. It should be noted that the entire closed-loop pattern of the multi-sinusoid 
io spectral observer is identical to the known in the art closed-loop resonator bank 
(CLRB) configuration, that is known per se (see, for example, N. Chenioguz, 
"'Tracking and recovering multiple sinusoids by an adaptive resonator bank" Proc. 
Fourth IASTED International Conference Signal and Image Processing, 2002, PR 
49-54). 

15 Fig. 2 illustrates a block-diagram of the CLRB configuration. In the 

mathematical model corresponding to the CLRB configuration, the /-th resonator Ri 
can be represented by 

In turn, the transfer function y-*yi , connecting the total input with the 
20 particular /-th output^ (or, equivalently, c,) can be obtained by 



(7) 



It should be appreciated that when the frequency of the transfer function 
(TF) achieves the location of the corresponding resonant peak, then the other 
resonators in the denominator may be ignored, and the magnitude of the TF for 
25 the corresponding sinusoid tends to 1 . 

It s hould b e r elevant t o n ote h ere t hat i n t he single tone case (i.e., when 
jV-1), the above multi-tone decomposition scheme reduces to a trivial zero-phase 



bandpass filter (BPF), where the corresponding transfer function Wi can be 
obtained by 

W=r ^ — (8) 

1-2(1 -0.5/)C /2 H +(l- r )z- 2 

It should be also noted that since the iV-sinusoid spectral observer 

represented by the system of equations (3) describes the propagation of both the 

in-phase and quadrature components, then the desired amplitude and phase of all 

the N sinusoids can be calculated. In particular, the corresponding quadrature TF of 

the z-th sinusoid can be obtained by 

Moreover, by utilizing the in-phase and quadrature components c/ and s/, 
respectively, the amplitude and phase of corresponding /-th sinusoid can be 
obtained by 

.A^JefTsf.. (10a) 

and 

q>. =atun(c t ./s f ). (10b) 



In order to take into account the fact that signal generated by the 
photodetector can be corrupted by optical disturbances and/or signal trends, the 
auxiliary states have to be included into the mathematical model described above 
for each disturbance and/or trend. Thus, a single-component bias term x^s has to 
be added to the mathematical model described above for the trend. In such a case, 
the one-state equation 

*«» (» + !> = T aUl9 x Bies (n) + G BJas e{n) (1 la) 

should be added to the W-sinusoid spectral observer represented by the system of 
equations (3) with the identity transition function TbiwH and gain Gnias. Here, the 
filter gain can be factored by a relaxation parameter Gb^s appropriate for the 
trend change. The observation function h can be updated by the component /*Bias = l- 
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Lilccwisc, a Nyquist single-state term x^si corresponding to the Nyquist 
noise component can be taken into account. In this case, the one-state equation 

x nm ( n + ') = T w,*^ (») + Gh4*<*) ( 1 lb ) 

should be added to the iV-sinusoid spectral observer represented by the system of 
5 equations (3) with the transition function 7)^. vl = -1, and the gain G^ m . The Nyquist 
gain has a negative sign and can be factored by a relaxation parameter y., so that 
Cijkts- -y. The observation function h can be updated by the component h^t m 1 . 

It should be noted that in the presence of periodic noise components, the 
//-sinusoid spectral observer should be updated by the proper resonators in the form 
i o of Equation ( 1 ) based on the knowledge of the noise component spectra. 

In the case when the central resonator frequency is unlaiown, the stationary 
CLRB version described above can be updated, in order to enable the estimation of 
the frequency-related parameter- For example, the adjustment of the frequency can 
be carried out by computing the phase increment over the one-step update, namely 
15 o>, (a) = pco } (// - 1) + (l - p)y> } (n) - f Pi (n - 1)] (12) 

where (pin) is the phase computed in accordance with Eq. (10b), and properly 
unwrapped, when required; and p is a relaxation factor which is a tuning parameter. 
The unwrapping of the phase assumes complementing the phase, when it jumps 
higher than 71 radian, so that the phase argument would maintain its continuity. The 
20 frequency adjustment translates the CLRB configuration into an adaptive resonator 
bank (ARB) form. 

Referring to Fig* 4, main units of the system of the present invention is 
illustrated, according to another embodiment of the invention. A ccording t o this 
embodiment, a synchronization is provided between the illumination assembly 131 
25 and the adaptive resonator bank (ARB) 143 operating in one of the regimes 
selected from a short serial mo3e, short serial-parallel mode and mixed-rate 
short-serial mode regimes. 

In operation, each light source (hereinafter all types of the possible light 
sources for conciseness will be referred to as LEDs) of the illumination assembly 
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131 turns on and off at a presumed rate producing periodic, nearly square probe 
optical waveforms. The total light signal (transmitted, scattered or reflected) is 
sensed by the photodetector 132, The multiple channel signal generated by the 
pholodetector 132 is fed to the analyzer where it (after conditioning and 
5 amplification) i s s ampled, e .g., at the rate o f about 2 50-300 kHz, and thereafter 
decimated. According to the invention, the signal after the sampling and decimating 
is treated by the DSP 141 in a multi-stage signal processing representing a 
particular demodulation-decimation technique. Examples of the demodulation- 
decimation will be described below in detail. 
10 As describe above, according to one embodiment of the invention, LEDs 

can generally operate in a short serial-parallel mode regime. In this regime, all the 
LEDs arc divided into several parallel groups, each associated with a particular 
modulation rate. In each group LEDs are involved in a short serial mode, one by 
one, for a short interval comprising one (or few) cycles of ignition. The number of 
15 the parallel groups, rates and intervals of the LED activity are design parameters, 
which should be adjusted in accordance with the system requirements. 

Alternatively, according to another embodiment of the invention, LEDs can 
generally operate in a mixed-rate short serial mode regime. In this regime, two or 
more generator-induced sequences with different modulation rates are initially 
20 mixed, and then the aggregated pulse train is fed to each LED in the short serial 
mode regime. 

According to one embodiment of the invention, the ARB can reserve one 
particular resonator for each corresponding LED (i.e., wavelength)- In other words, 
only fundamental harmonics are extracted from the composite signal. Thus, for 

25 given N LEDs (i.e., N optic channels), N corresponding resonators are to be 
involved, according to this embodiment. In such a case, the resonant frequency of - 
each resonator fits the modulation rate of the corresponding LED. It should be 
noted that when the LEDs are activated concurrently, in a parallel mode, their 
frequencies should be different, while when the LEDs are activated serially, the 

30 modulation rate for all LEDs may be identical. 
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It should be appreciated that when required, several resonators (R^ through 
R mm ) can be used for each optic channel, corresponding to different harmonics 
associated with this channel (see Fig* 3), 

As shown in Fig. 3, the synchronization between the MUX and the ARB is 
5 provided. For example, when a particular LED is turned on, the input composite 
signal is fed to the appropriate resonator, and the resonator's output is linked to the 
total output signal. In turn, when a particular LED is turned down, the 
corresponding resonator's input is locked and the resonator's output does not 
contribute into the total output. However, the resonator with the locked input can 
io sustain oscillations in the autonomous mode, in order to be synchronized in time 
with all the input components. Hence, the resonator can produce a contemporary 
signal, irrespective of the LED's activity. 

According to the invention, the adaptive resonator bank can be further 
designed to handle slowly varying components of the composite signal provided by 
15 the photodctcctor 132, e.g., signal trends. The signal trends are usually described by 
its first derivatives, e.g., the dc (constant component), slope (rate). Optionally, the 
signal trends can be described by higher signal derivatives, e.g., acceleration and 
jerk. A corresponding portion of the ARB responsible for handling the trends is 
referred to as a trend filter. According to the invention, the trend filter is attached to 
20 the CLRB configuration (see Eq. (11a)) with the corresponding trend filter F Tr md 
for the_signal trend. 

It should be noted that the number of the derivatives used in the trend filter 
depends on the trend's order. For instance, because in practice the dc is the most 
dominant component, the appropriate trend filter has the first order. In turn, the 
25 trend filter, relying on the dc component and the slope, can have the second order, 
etc. Hence, the trend filter has a predetermined order and M capable to track the 
signal trends, irrespectively of the LED's modulation mode. 

Likewise, a further portion of the resonators attached to the CLRB 
configuration can be responsible for the rejection of different periodic disturbances 
30 having either known or unknown frequencies. Such resonators can have 
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predetermined in advance resonant frequencies, when the interference location in 
the frequency domain is known. When the frequencies are unknown, the 
frequencies can be properly adjusted, for example by using the scheme described 
by Eq. (12), Thus, the application of the CLRB configuration (represented by the 
5 spectral observer scheme) enables to eliminate various narrow-band disturbances 
from the input composite signal, and therefore prevents their influence on the 
LED-related components. In the same manner, resonance filters., related to other 
types of interference, may also be introduced into the ARB framework. 

It should be relevant to emphasize here that instead of eliminating various 
10 narrow-band disturbances from the composite input signal by applying a notch 
filter before deriving LED's related signals, the present invention teaches to 
eliminate the interference spectral components simultaneously with the recovering 
of desired components. This processing scheme is common for the LED's related 
signals as well as for the trends, disturbances and noise. According to the invention, 
15 the elimination of the trends, disturbances and noise is carried out by applying a 
corresponding ARB. Such an approach provides an accurate separation of the 
desired channels from the undesired spectral components, whereas a "bandstop" 
rejection of the undesired disturbances by the notch filtering before the 
decomposition of the LED's induced signals may affect the desired spectral 
20 components as well. 

For example, the signal processing by the ARB shown in Fig. 3 can be 
formulated as follows. The treated signal y(t) can be represented as a superposition 
of induced by LEDs' periodic multiple-harmonic signals, slow trend signal, 
sinusoidal or periodic disturbances, and a noise. 
25 In a general scenario, the LEDs can be driven by the r pulse trains having 

different pulse rates co r , /— 1, 2, ... The /-th LED is activated by the -corresponding 
/-th portion of the pulse train (if the train sequence is single) or by a superposition 
of the /-th portion in each r-th train (if several triggering train sequences are applied 
at a moment of time). In such a case, the corresponding composite signal provided 
30 by the photodetector can be presented by: 
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N ^ 

l-r\ r =1,2....* =1.3...- (13) 

+ J D. sin(a>/ + ^ )+ Trend (t) + v(/) 
/=j 

The first sum-term in the right-hand side of Eq. (13) represents a sum of the 
LED-rclatcd waveforms over (i) N LEDs, (ii) fundamental to higher harmonics 
associated with particular pulse sequence waveform, and (Hi) all the driving pulse 
5 sequences involved. The second sum-term relates to the sum of sinusoidal 
disturbances. The next term, denoted as Trend, refers to a slow time-varying trend, 
while the last term designates the white noise. 

The link between the LED activity and resonators is established via the LED 
switch function %J. The factor %J ^xfft) * s the binary (0/1) function that equals 1, 
10 when the /-th LED is activated by the r-ih pulse train; or 0, otherwise. In this 
notation, the superscript indicates the /-th pulse position in the pulse train driving 
the /-th LED. As the LEDs are driven by a s ingle pulse sequence, i.e,, r=l, the 
subscript r in the LED switch function may be omitted, resulting in a reduced 
notation;}/ =^ '(/), w here / denotes the LED, M,.. JV. As the number of pulse 
15 sequences involved is more than 1, the LED switch function obtains a subscript r y 
2,... denoting the driving sequence. 
For treating the signal described by Eq. (13), for example, the Kalman-Hke 
spectral observer (SO) scheme, described above, with "nearly" optimal gain can be 
employed. Such a modified SO scheme takes into account that the particular /?-th 
20 resonator for a single tone cor can be involved only at particular instants of the 
LEDs 7 activity indicated by the binary (0/1) switching function Xc > to wit: 
*,(* + !) = r,x,(n) + X l(n)G,e(ii) (14) 
In other words, the prediction step in Eg. (14) is completed as usual (see,,- 
Eq.(l)), while the update step is controlled by the LED's switching function %} . 
25 Likewise, the spectral observer is a generalization of the single-tone version, 

to wit: 
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x(n + 1) = 7*(w) + %{n)Ge(n) 9 ( 15) 
e{n) = y{n)-%{n)hx{ri) 

where ^(n) should be properly composed from the binary switching functions 

Xr( n X ^ ,2,...; /~1 , An example of the detailed description of the switching 

function i(n) will be given below. It is relevant to note here that the computation 

5 the prediction error also depends on the function x(n) and should be in agreement 

with the resonators states. Further, it is found convenient to aggregate the functions 

%(ri) and h into a resulting observation function h(n) composed of the particular 

LED's switching functions yj(n). 

It should be appreciated by a person versed in the art that the state vector x 

io should, minimally, comprise the states associated with the LEDs 5 in-phase and 

quadrature components. In a more general consideration, when periodic 

interference (optical and/or electromagnetic) disturbs the signal, the corresponding 

auxiliary states should be involved into the model. 

For example, an electrical interference at a fixed frequency (e.g., 100 Hz) 

j 5 may affect the PD signal. This continuous-time interference may appear 

irrespective of the LEDs' control mode. Hence, the resonator operating at the 

corresponding frequency of 100 Hz should be included into the SO scheme with 

corresponding switching function, that is always equal to L 

Moreover, when the value of the frequency is unknown (or slightly varies 

20 with time), a frequency estimator can be applied, e.g., based on Eq. (12). For this 

puipose, the corresponding frequency estimator can be updated recursively, as 

indicated above. 

Further, three examples of the configuration system of the present invention 
will be illustrated, 

25 

EXAMPLE 1 

Fig. 4 illustrates a schematic diagram of a system operating in a short serial 
mode regime. According to the invention, the MUX 12 is configured to operate in a 
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short serial mode regime. Fig. 5 illustrates an example of waveforms produced by 
different elements of the system shown in Fig. 4. 

According to this example, a series of activating pulses F (provided by the 
generator 11) are fed to the MUX 12- The MUX 12 operates in such a manner that 
5 only one on-off ignition cycle is provided for the modulation of each light source 
131. The pulses F 1 , F 2 , ...and F N correspond to the operation regime of the light 
sources having numbers 1, 2, ... and N, respectively. Hence, the light sources are 
run in series, one after another, each LED over one ignition circle. Further, the 
pulses indicated by symbol PB correspond to the composite signal generated by the 
10 photodetector 132 in response to the light provided by the light sources 131 
activated by the pulses F 1 , F 2 , ...and F N . The analog signal PD, converted by the 
ADC 140 into a high-sample-rate digital signal, and re-sampled by the first 
declinator 142 to a lower rate SRi, is then decomposed by the ARB 1 43 which 
produces a set of pulses associated with resonators R 1 , R 2 , ...and R N , 
is corresponding to the light sources having the numbers 1, 2, N 9 respectively. 
Finally, the pulses can be filtered out and the output signals 1,2,..., jV representing 
light intensity E corresponding to the light sources having the numbers 1, 2, . . N 9 
respectively, can be obtained and used for determination of blood parameters. 

For instance, when a number of the light sources equals 10, a sampling rate 
20 of 250 kHz can be used for sampling the signal PD generated by the photodetector. 
Preferably, the modulation frequency for the activation of the light sources (e.g., 
LEDs) is synchronized with the operation of the ARB. Thus, the light modulation 
frequency can be selected within a wide range restricted from the high frequency 
limit by the system's time constant, and from the low frequency limit by the low 
25 frequency disturbances and their most essential harmonics. In the short serial 
. regime of this above example, each LED can be turned on over 40 samples, and 
turned off over next 40 samples. Hence, the LED's ignition period is 80 samples, 
which is equivalent to the modulation rate 250/80-3.125 kHz. 

It should be noted that although triggering the light sources provides nearly 
30 ideal square pulses (see the pulses F 1 , F 2 , ... and F N ), the signal produced by the 



photodetcctor (that is usually integrated with an amplifier) represents a periodic 
signal with a rather complex shape and, accordingly, with a wide harmonic content. 
Thus, a decimation slep should avoid the aliasing phenomenon, i.e., a leakage of 
the higher-order harmonics into the lower frequencies. 
5 It should be relevant to note here that the length of conventional 

anti-aliasing low-pass filters is usually several (e.g., 3-4) times higher than the 
decimation ratio and this may cause undesired coupling between the LEDs. 
Therefore, the decimation performed with the help of a conventional anti-aliasing 
filter is not appropriate here, due to the strong necessity to keep independence 
10 between the different LEDs. Nevertheless, the method of the present invention 
allows to reduce the 250 kHz sampling rate by avoiding the aforementioned 
straightforward decimation procedure. The reduction can be done by relying on the 
particular spectral content of the given signal. 

More specifically, noting that the generated signal is nearly an ideal 
15 symmetric waveform, the odd. harmonics of the spectrum can only be taken into 
account, while the even harmonics, describing the pulse asymmetry, can be ignored. 
By integrating consequent groups of the samples, it is possible to reduce the 
sampling rate and narrow the spectrum bandwidth, by keeping in mind the location 
of the desired odd and undesired even harmonics in the frequency domain. 
20 In particular, one should take into account that the high-order odd harmonics 

can merge, after the decimation, with the lower odd harmonics, while high-order 
even harmonics can merge with the lower even harmonics. Thus, by using the 
two-sample summation, the sampling rate can be reduced to 125 kHz. In such a 
case, the even 80-th harmonic merges with the dc, the 79-th harmonic merges with 
25 the 1-st (fundamental) harmonic, etc. Therefore, after the two-sample summation 
stage, the only AO harmonics remain. Likewise, the summation of several samples, 
may result in the situation when only few lower harmonics remain in the 
compressed spectrum. Thus, only small number of resonators can be required for 
each channel For instance, the summation by 10 samples results in a spectrum 
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comprising only 1-st and 3-rd harmonics. In other words, only two resonators are 
required for each channel. 

Accordingly, the s ummation by 20 samples results in such a situation, in 
which all the odd harmonics would merge with the 1-st harmonic at 3.125 kHz. In 
5 turn, the even harmonics would partially merge with the dc component, and 
partially with the Nyquist frequency, e.g., 6.25 kHz. This approach provides a 
spectrum with a single odd harmonic, which needs only a single resonator. 

Thus, the summation by 20 samples reduces the original sampling rate of 
250 kHz to the rate SRi of 12.5 kHz (see Fig. 5). Accordingly, the integrated 
io sequence of pulses contains a 2-sample on-portion and a 2-satnple off-portion. 
Hence, the total sample rate of the sequence is set to 3.125 kHz, i.e., 4 samples per 
one on-off ignition cycle. Accordingly, the basic frequency of 3.125 Hz can merge 
with all the odd harmonics presented in the original (6.125 kHz width) spectrum, 
i.e., with the harmonics corresponding to 3.125x(3, 5,..., 39) kHz. In turn, the 40-th 
15 harmonic falls onto the dc component, the 38-st - onto 6.25 kHz (Nyquist), the 
36-st - onto the dc, the 34-st - onto the Nyquist signal, etc. As noted above, only 
the odd harmonics are taken into account, while the even harmonics are ignored. 
Hence, the resonator should match the 3.125 kHz, while the dc and Nyquist 
components may be rejected. 
20 According to the invention, the MUX distributes the pulse sequence 

between the different LEDs. Accordingly, the 1-st LED receives the first pulse of 
the sequence, the 2-nd LED receives the second pulse, etc. (see, for example, the 
pulses F*, /=1, 2,...,Ar shown in Fig. 5). Hence, each LED, in turn, can be switched 
on and off over a period of 4 samples (at the 1 2.5 kHz rate), and. then would remain 
25 passive over the next 9x4=32 samples. The minimal shift between the consequent 
active LEDs is 4 samples. - 

The photodetector integrates the optical signals induced by all N LEDs and 
generates the composite signal (PD in Fig. 5). As can be seen from the shape of the 
PD signal, the square-waveform of the original LED's pulses is slightly 
30 disturbed by the PD inertia. It should be noted that the particular pulse frequency of 



-26- 

3,125 kHz can be selected to provide a minimal leakage between the consequent 
LEDs' waveforms. 

The signal PD is fed to the ARB for filtering and channel separation. The 
ARB corresponding to this example can comprise at least N (e.g. N=1Q> as in the 
5 present example) resonators having identical fundamental frequency 3.125 kHz, 
and producing R 1 , R\ R N signals. All these resonators operate in parallel, 
however at each moment the PD signal is fed to only one of them, owing to the 
LED switch functions % . % , which are binary 1 or 0 indicators synchronizing 
the ARB with the LEDs 7 operation. Hence, for the given 10 wavelengths, the state 
10 vector of the spectral observer comprises 20 states (one in-phase and one 
quadrature for each resonator), where two parameters in the state-space vector are 
reserved for the dc component (0 Hz) and the Nyquist frequency (6.25 kHz), 
respectively. 

Each output of resonators represents a periodic signal having a certain 
15 amplitude. The shape of the signals R 1 , R 2 , . . R* can be different from the square 
waveform of the original LED's pulses F*. It should be appreciated that when 
required, the original pulse form can be reconstructed with the help of auxiliary 
resonators related to higher harmonics. 

At the last stage, when required, the 12.5 kHz sampling rate can be further 
20 reduced by using a conventional decimation, e.g., to the 25 Hz sampling rate, where 
ail the channels would indicate the time evolutions of the light responses 
synchronously. 

EXAMPLE 2 

25 Referring to Fig. 6 and Fig. 7 together, an example of the system operating 

in a short serial-parallel regime .and the corresponding waveforms are illustrated. 
The system of this design is different from that described in Example 1 in that the 
ten LEDs are divided into two parallel groups. The number of the LEDs in each 
group may be arbitrary, e.g., 6 and 4, or 9 and 1, or 5 and 5, etc. The first group 

30 comprises the LEDs %\ through Xi M 5 and the second group comprises the LEDs \z 
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through Xi \ where M±L=*10» The first and second groups are fed by MUX1 and 
MUX2, respectively, and operate in the short serial mode each, however at different 
rates. For example, the rate of 3.125 kHz for the first group and the rate of 2.5 kHz 
for the second can be selected. 
5 Accordingly, the Generator 11 outputs two pulse sequences Fi and F2 shown 

in Fig, 7. Contrary to Example 1, in the present example the 250 kHz samples are 
integrated by groups of 10 samples, thus reducing the sample rate from 250 kHz to 
25 kHz. Thus, in the first sequence each activating pulse contains 8 samples per 
on-off cycle (4 samples for the on-portion and 4 samples for the off-portion), while 
10 each pulse in the second sequence contains 10 samples per on-off cycle (5 samples 
for the on-portion and 5 samples for the off-portion). It should be noted that owing 
to the sampling rate of 25 kHz, the signal spectrum bandwidth is 12,5 kHz, i.e. half 
of the sampling rate. 

To this end, the spectrum having the width of 1 2.5 kHz can exhibit only the 
15 1-st and 3-rd (9.375 kHz) harmonics of the first, i.e., the 3.125 kHz rate sequence. 
With regard to the second, i.e., the 2.5 kHz rate sequence, the situation is similar. 
The 5-th harmonic of the fundamental frequency of 2.5 kHz (i.e., 12.5 kHz) can 
merge with the Nyquist frequency and should be disregarded. Therefore, only the 
1-st and 3-rd harmonics (7.5 kHz) of the fundamental frequency of 2.5 kHz can be 
20 used. Note that using the 3-rd harmonic can increase the system bandwidth, when 
compared to the previous example. Thus, two resonators should be reserved for 
each LED in the corresponding ARB. In this case, the total LED output is obtained 
by aggregating the 1-st and 3-rd harmonics. 

Accordingly, the pulses from the first sequence Fi are fed to the 
25 corresponding LEDs in the fiist group. In particular, the subsequence Fi 1 drives the 
LED A,t\ the subsequence Fi 2 drives the LED fa 2 , etc. In turn, the pulses of the' 
second sequence F 2 are fed to the second group of the LEDs. In particular, the 
subsequence F2 1 drives the LED Xz\ subsequence F 2 2 drives the LED Xi 2 , etc. It is 
relevant to recall that in the notations used here the subscript indicates the 
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Generator sequence (in particular, the first or second), while the superscript - the 
pulse in the sequence (1 to N). 

The signal PD in Fig. 7 represents an output waveform generated by the 
photodetector in response to the signals induced by the LEDs. The signal PD is 
5 then decomposed by the ARB 143. According to this example, the resonators of the 
ARB operate in a mixed time-frequency separation mode. Thus, the resonators 
associated with the 1-st group of the LEDs operate in series, A part of these 
resonators is synchronous with the basic frequency 3.125 kHz, and a remaining part 
is synchronous with its 3-rd harmonic having the frequency of 9.375 kHz. 
to More specifically, each LED is associated with two resonators or, in other 

words, with two pairs of states in the spectral observer state vector, wherein a first 
pair (in-phasc and quadarture) corresponds to the 1-st harmonic and a second pair 
corresponds to the 3-rd harmonic. The spectral observer state vector comprises, 
therefore, 40 states for 10 wavelengths. 
15 For instance, a two-component state vector x% can relate to the in-phase and 

quadrature of the 1-st harmonic of the 1-st LED. The corresponding resonator is 
specified by (i) the appropriate LED switch function %i l (related to the 1-st 
sequence and 1-st pulse position in the 1-st sequence) and (ii) the appropriate gain 
and transition matrix (matching the rate 3,125 kHz). Accordingly, a two-component 
20 state vector j can relate to the in-phase and quadrature of the 3-rd harmonic of 
the 1-st LED and the same switch function while the gain and transition matrix 
should be in agreement with the tripled modulation rate 9.375 kHz. Likewise, the 
next two states (four components) are associated with the 1-st and 3-rd harmonics 
of the 2-nd LED synchronized by the switch function and the same for all M 
25 LEDs of the first group of LEDs. 

For the L LEDs of the second group, the resonatprs^should be in agreement 
with the modulation rate of 2.5 kHz and its 3-rd harmonic. It should be appreciated 
that the corresponding modification can also affect the elements of the SO 
transition matrix T and the gain vector G. The resonators of the second group are 
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synchronized with second sequence of pulses through the LED switch f unctions 

Another difference of Example 2 from Example 1 is in that the output of the 
resonators i s obtained by aggregating the 1 -st and 3 -rd harmonics related to the 
5 same LED, Accordingly, the algorithm should contain a step representing the 
aggregation of the resonators corresponding to different harmonics of the same 
LED signal. 

Turning back to Fig. 7, examples of the outputs of the ARB obtained by the 
aggregation of the 1-st and the 3-rd harmonics are depicted. In particular, Rj 1 is 
10 associated with the LED A,i\ an output R2 1 is associated with the LED Aa 1 , etc. 

It should be noted that when the mixed serial-parallel mode regime is 
employed, the LEDs' silence intervals become smaller than in the pure serial mode 
regime of Example 1 , thus the smoothing ability and the bandwidth of the system 
are improved. Moreover, the modulation rates 2.5 Hz and 3.125 kHz are set 
15 sufficiently far from each other, thus preventing the coupling between the LEDs. 

The number of the LEDs in the parallel groups may be arbitrary selected to 
meet the system requirements. For example, one particular LED, providing a 
reference wavelength, may be operated continuously at the modulation rate of 
2.5kHz, while the other 9 LEDs may operate in a short serial mode. 

20 

EXAMPLE 3 

Referring to Fig. 8 amd Fig- 9 together, another example of the system 
operating in a mixed-rate short serial mode regime and the corresponding 
waveforms are illustrated. According to this example, the modulation mode regime 
25 is different from those used m Examples 1 and 2. 

The Generator 11 provides two pulse sequences Fj and F 2 at-rates of 3.125 
kHz and 2-5 kHz, respectively. The M UX1 and MUX2 subdivide each of these two 
sequences imo ten subsequences. In particular, the sequence Fi corresponds to the 
signals Fi\ Fj 2 ,.„, Fi N , while the sequence F2 corresponds to the signals F 2 \ 
30 F2 2 ,, . F2 N . 
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The first short serial mode regime (the pulses Fi\ Fi 2 ,..., F^) is applied to 
the N LEDs (where JV^IO) at the modulation frequency of 3.125 kHz, as in 
Example L The LED switch functions associated with the first sequence are 
denoted as %\ Z^- At the same time, another short serial mode regime (the 
5 pulses F 2 \ F 2 2 ,. . ¥ 2 ") is applied to the same N LEDs having the modulation rate 
of 2.5 kHz. Accordingly, the LED switch functions associated with the second 
sequence is denoted as / yz N - In other words, the superposition of the pulses 

Fi 1 and F 2 * drives the LED Xt, the superposition of Fj 2 and F 2 2 drives the LED %2> 
etc. Finally, the superposition of the pulses Fj 10 and F 2 10 drives the LED A,io. 
10 According to this scenario, each LED is activated by a combination of two 

modulation sequences. This case is different from that of Example 1, because each 
LED receives a combined complex-form pulse sequence, rather than a simple 
square-pulse sequence. 

At the receiver side, for the decomposition of the PD signal, the number of 
15 resonators properly increases. In particular, each recovered signal R/ (r=l,2; 
/=1,..., N) needs two resonators, for the 1-st and the 3-rd harmonics, respectively. 
Each LED is triggered by the two subsequences, and therefore each channel needs 
four resonators (pairs of states), such as a resonator for the 1-st harmonic 
corresponding to the modulation frequency of 3.125 kHz, a resonator for the 3-rd 
20 harmonic of 3.125 kHz, a resonator for the 1-st harmonic of 2.5 kHz, and a 
resonator for the 3-rd harmonic of 2.5 kHz. Aggregation of these four resonators 
provides a total output signal for the corresponding channel. 

Each resonator operates in the manner described above, under the control of 
the LED's switch function. The output signal is aggregated from the in-phase and 
25 quadrature of each resonator (as in Example 1), then it aggregates the I -st and 3-rd 
harmonics (as in Example 2), and, finally, it combines outputs associated withHhe 
modulation pulses at different rates of 3.125 kHz and 2.5 kHz, correspondingly. 

It should be appreciated by a person versed in the art that driving each LED 
by superposition of two modulation sequences can increase the measurement 
30 quality. This improvement can be achieved exclusively at the expense of moderate 
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increase of computations, whereas the hardware remains the same. It is noteworthy 
that the mixed-rate scenario of Example 3 allows a reliable and accurate extraction 
of the signal components without significant interaction between the two 
concurrent frequencies of 2.5 kHz and 3.125 kHz. 
5 It should be appreciated that despite the differences of the above examples, 

all of them can be implemented within the same DSP framework. A generalized 
algorithm implementing the method of the present invention complying with these 
three examples and other scenarios for decomposition of multi-channel signals is 
presented as follows. 



Design Parameters: 

Slow Convergence Parameter, £«, gy 

Fast Convergence Parameter^ y (U yj 

Frequency Convergence Parameter, p a , pj 
1 5 Forgetting Factor ju, 

Fixed- Frequency Disturbances d=\,2^..\ 

Unknown Disturbances Initial Guess mo, u=l, 2,. . .; 

Modulation Rates cor> r = 1 , 2, , . . ; 

Indices of LED s /=1, 2,..., N\ 
20 Indices of Harmonics #=1,3,...; 

LED Switch Functions %*, 2,.-.; 7-1, 2,., .,7V; 

State VeclorDesignationi 
Xi - dc (single state) 
25 X2 - Nyquist (single state) 

x 3 - l" Modulation Rate, I*' LED, I st Harmonic (in-phase&quadrature states ) 

x; - r Modulation Rate, f LED, k' h Harmonic (in-phase&quadrature states) 

30 x P - a' known-frequency disturbance m (in-phase&quadratue states) 

- u lh unknown-frequency disturbance qju (in-phase&quadrature states) 

35 Transition Matrix of i" Resonator for t H LED. r th Modulation Rate & k' h 
Harmonic 

C,=cos(kG)r), $i*=sin(kcor), 
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C, S, 

-s, c, 

c, s, 

i=ir.u. 2...; k=L 3,...; 1=1,. ...N 



Transition Matrix ofp' f> Resonator related to d h jvxed-frequency Disturbance 
C p =cos((oj), S p =sin(aya), 



10 



15 



20 



25 



30 



35 



p=pd. d-l, 2.... 

Fixed Part of Integrated Transition Matrix 
T 0 =\@(-l)6bd(7,. Mr*, r-1. 2...; /c=l. 3,. 
2,...) 

Fixed Part of Gain 

i"/r P ii/, r ^\> 2..-; 3,..,; l=*l,„.,N 
p=p<t> d=l,2 9 ..r K 



Initialization for n-O: 

Fast Relaxation Factor y(0)=y d 

Slow Relaxation Factor e{0)-€a 

Frequency Convergence Factor f{0)- p a 

Unknown Disturbance Frequency: cou(0)=g> U u. 2,...; 

Compute For Time Instant n=l. 2. . . . 

Current Observation Function 

A=A(/z)=[l, -1, 0, .... ^ 0,...,1, 0.....1, 0] 

Current State Update Function 



l=\,...,N)®t>d(T p9 p=p dt d=l, 



Transition Matrix of Resonator related to u h Varying-frequency Disturbance 
Q= cosf&J, S q = sin(ou), 



To= 



40 w-1, 2,.,. 



Integrated Transition Matrix (Updated) 
T=T 0 @bd(T q , q=cf u , u=l, 2,...) 

5 Updated Gain 

G- [Go, s( C q , S< h ...)], 4=4* u=\ f 2,...; 

x(n+\\\\)=Tx(n) - Prediction 
e(n)=y(n)-h(n)x(n+ 1 \n) - Prediction Error 
10 x(n+\)=x(n+\\n)i- %(n)Ge(n) - State Update 

Amplitude Related to f LED, /=], ... f N 

(Computed From Tn-Phase/Quadrature Components and Used Harmonics) 
15 Ai=$qrt(I£f+&j* f i=i,. k j, r=»l, 2...; k=l, 3,...; 

Unknown Disturbance Phase 
(p f/ (n+])=arctan(c</set)+ k q n, q=q u , «=1, 2,...; 

£ v w aw unwrapping factor preventing higher than n jumps in <pq(n+\) 
20 compared to the preceding term (pq(n). 

Unknown Frequency Update 

#//z+i;=p^(n)+(l-^ q=q„, u=l, 2,..,; 

25 Relaxation Factors: 

y^jtff+Cl-tth/h e=>{&+(l-/4)ei, pr=}ip+(\-ij)p\ 

End of the loop 

Tn the above algorithm, the symbol © denotes block diagonal matrix 
30 expansion operator, while the symbol bd(*) relates to the block diagonal matrix 
composed of its parameters. 

As such, those skilled in the art to which the present invention pertains, can 
appreciate that while the present invention has been described in terms of several 
embodiments, the concept upon which this disclosure is based may readily be 
35 utilized as a basis for the designing of other structures, systems and processes for 
carrying out the several purposes of the present invention. 

It will also be understood that the system according to the invention may be 
a suitably programmed computer. Likewise, the invention contemplates a computer 
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program being readable by a computer for executing the method of the invention. 
The invention further contemplates a machine-readable memory tangibly 
embodying a program of instructions executable by the machine for executing the 
method of the invention. 
5 Also, it is to be understood that the phraseology and terminology 

employed herein are for the purpose of description and should not be regarded as 
limiting. 

in the method claims that follow, alphabetic characters and Roman numerals 
used to designate claim steps are provided for convenience only and do not imply 
10 any particular order of performing the steps. 

Finally, it should be noted that the word "comprising" as used throughout 
the appended claims is to be interpreted to mean "including but not limited to". 

It is important, therefore, that the scope of the invention is not construed as 
being limited by the illustrative embodiments set forth herein. Other variations are 
15 possible within the scope of the present invention as defined in the appended 
claims and their equivalents. 



